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I. INTRODUCTION 



The problem of back-reaction arsies in a number of contexts ranging from 
the modehng of the pre-equihbrium phase of quark-gluon plasma production [1,2] 
to studies of inflationary cosmology [3-5]. In general terms, back-reaction consists 
in a situation where a system can be viewed as governed by the mutual interaction 
between a field, usually taken as classical, and pairs of charged particles that are 
produced through its presence. In one version, the initial condition posits the 
pre-existence of the field and an approximate vacuum insofar as the particles are 
concerned. Pairs can then tunnel out of this vacuum by reason of their interaction 
with the field in a purely quantal process. These pairs are accelerated by the field, 
producing a current which in turn acts back on the field. The back- reaction leads 
to plasma oscillations as one of its key characteristic effects. 

Within the past few years, it has proved possible to carry out detailed 
calculations for back-reaction in which boson [6] or fermion [7] pairs are produced, 
including the case where boost-invariant coordinates are used [8] as is appropriate 
for the study of the quark-gluon plasma. (See also the didactic article in Ref. [9] 
and the review in Ref. [10].) Somewhat surprisingly, it has emerged that calcula- 
tions [1,2] based on the much simpler transport formalism follow the field-theory 
results in detail provided that quantum fiuctuations are smoothed and quantum 
statistical effects are introduced into the transport equations [6-10]. That is, it 
proves possible to mimic the quantum tunneling through the use of a source term 
based on the Schwinger mechanism [11] in a transport equation. The squares 
of the quantal mode amplitudes, describing the momentum distribution of the 
produced pairs, even agree remarkably closely with the corresponding transport- 
function momentum dependence once smoothing and statistics are included. There 
are also formal derivations [12,13] of the mapping from field theory to transport 
formalism in this context, providing [13], for example, an understanding of how 
nonmarkovian features enter the theory. Thus the back-reaction problem provides 
a model of a system sufficiently simple that the nonperturbative field-theory de- 
scription can be fully solved numerically and approximations to it — notably the 
transport formalism — can be studied in detail, but with enough complexity to al- 
low for a richness of phenomena and features. This is especially interesting when 
the nonequilibrium thermodynamics of back-reaction is considered [8] . The field- 
theory version then refers to the time evolution of a pure quantum state, while 
the transport formalism has a well-defined entropy which increases with time. 
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One may also hope that essential features of the quark-gluon plasma lack- 
ing in this simple back-reaction model will eventually be incorporated and their 
effects studied. These must surely include (i) the nonabelian nature of quantum 
chromodynamics which governs the plasma; (ii) the nonclassical nature of the 
chromoelectric field that enters there; and (iii) the fact that the plasma lives in a 
finite volume, not in the infinite volume treated so far in back-reaction. To the 
degree that a mapping from field theory to transport formalism can be found, say, 
in the presence of each of these features separately, the use of the relatively simple 
transport formalism for the description of the plasma gains considerable support. 
Towards this end, some progress has been made in the formulation of back-reaction 
for the nonabelian Yang-Mills case [14]; a detailed study [15] has been carried out 
for the lowest quantal correction to the classical electric field used so far; and it is 
the purpose of this paper to study back-reaction in a (transversely) finite volume 
of a geometry chosen for its eventual applicability to the quark-gluon plasma. 

We consider back-reaction in an infinitely-long cylindrical sleeve, as shown 
in Figure 1, for an abelian and classical electromagnetic field. The choice of a 
volume whose finiteness is only in the direction transverse to that of the electric 
field and particle current was made with a view to eventual extension to the quark- 
gluon plasma, where boost-invariant coordinates would be used for the longitudinal 
variables. (In one spatial dimension, the case of back-reaction on a finite line 
segment has been considered [16] in the past.) The symmetry here about the 
cylindrical axis is, of course, a sizable simplification for our problem, though it 
still leaves in it the complexities of mode mixing and of renormalization, and 
appears to be just barely amenable to handling at the numerical level with present 
100-MHz or 150-MHz computer processors. We develop below the formalism for 
back-reaction in this geometry, first for field theory and then for the transport 
formalism, then discuss briefly numerical methods and approximations for each, 
and compare results for the two approaches. 
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II. FIELD-THEORY FORMALISM 



We consider QED for scalar particles of charge e and mass m, governed 
by the Klein-Gordon equation 

[{d^' + ieA''){d^ + ieA^)+m^]^ = 0, (1) 

coupled with the classical Maxwell equation 

C>2 



5^2-V^J^^ = 0'")- (2) 

Here is the classical four-vector potential and j^^ is the four-current, taken in 
an expectation value for some initial configuration of the charges. The current is 
given by 

= lie{$t(aM + ieA^')^ - [{^^' + ^eA^')^]^^ 

2 (3) 

- + zeyl^)$]t + [{^^' + ieA>')^]<5^}, 

in terms of a second-quantized field $. 

We choose the geometry of an infinite cylinder of radius R (see Figure 
1) and restrict the initial conditions to be independent of the polar angle ip and 
assume that initially there are no regions of nonvanishing net charge j^{r,t) and 
that the initial field configuration consists of E(6, t) || z, where b = |b| is the 
radial coordinate. The current j(r, t) then remains parallel to the axis of the 
cylinder z. This is because (see Figure 1) the opposite charges of a produced 
pair are accelerated by E in opposite directions parallel to z. The magnetic field 
B(b, t) that is produced has solenoidal geometry, whence both the positive and 
the negative charges undergo transverse acceleration of the same direction and 
magnitude, canceling any transverse contribution to the current. Selecting a gauge 
such that A'^ — (0; 0, 0, A{b, t)) has a 2;-component only, the Klein-Gordon equation 
becomes 

r „ / ^ \ 2 1 

^b,<p,z;t) = 0, (4) 



where 



^ db'^ bdb b'^dif^' 



We choose boundary conditions such that there are no particles on or 
outside the cylinder wall at 6 = i?; the electric field is also taken to vanish at and 
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beyond h = R. These are boundary conditions appropriate to electromagnetism. 
One could of course choose bag boundary conditions in their place. In expanding 
$ over creation and annihilation operators, we then use the orthonormal basis set 



0«n(b) = — -7-7— ^-"Y^, (6) 
V2tt R Ji{k±[ln)R) 

where k±{ln) ~ zin/R, zin being the nth zero of the Bessel function of order I. In 
the following, we indicate the pair {Z, n} by n; the lack of (^-dependence restricts 
our interest to Z = in any event. Then 



where ak,n annihilates particles of longitudinal momentum k and radial mode n, 
b^-k n creates antiparticles with —k and n, and <^fc ^ = [k'^ + /c^(n) + m^]-*-/^. The 
usual commutators pertain, 

[afc,n(t), 4',n'(0] = [hAt), 4',n'(0] = 27r5(/c - k')6n,n' , (8a) 
KnW, 4',n'W]=0, (8b) 

and so forth. We must relate these operators to their values at t = 0, which we do 
by means of a Bogolyubov transformation [17] 

ak,n{^) = ^[Uk,nn'{t)ak,n'{Q) + Vk,nn' {t)b^_k,n'(^)i i^^) 



and 



where 



and 



bU,nit) - Y>l,nn'{t)^k,nm + 4,nn' (^)&^-fc,n' (0)] > (9b) 
n' 

«fe,nn'(0) = 5nn', ^^fe,nn' (0) = 0, (9c) 

J]K,nn(^)Wfe,n'n(^) " Vk,nn{t)vl^n'n{t)\ = ^nn' , (9d) 
n 

'^[uk,nn{t)Vk,n'n{t) - Vk,nn{t)Uk,n'n{t)] = 0, (9e) 



in an obvious extension of the usual [17] forms. 
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Defining the mode ampfitudes 

fk,nn'{t) = . 0^ .U2 ['^k,nn'{t) +Vl^^,{t)], (10a) 

with the initial values 

fk,nn'{Q) = JTr-Q—^^nn', /fc,nn' (0) = -i ( J Snn' , (10b) 

l^^/c,n^ \ ^ / 

the scalar field becomes 

^e^'^</>n(&)[/.,nn'(Oafc,n'(0)+/feV(^)^-fe,n'(0)]- (H) 

nn- 

The normalization of the mode amplitudes in Eqs. (10) is taken so as easily to 
accommodate the canonical commutation relations 

[$(r,t), ^(r^^)]=^5(r-rO, (12) 

where 11 = $ is the field conjugate to Equation (12), with the commutators of 
Eqs. (8), leads, among other similar relations, to the wronskian condition for the 
mode amplitudes, 

^[fk,nn{t)fk,n'nit) " fk,nnit)fk,n'n{t)] = iSnn' , (13) 
n 

which is guaranteed by the dynamics to be satisfied at all times if it is fulfilled 
initially, as it is of course by Eq. (10b). 

Substituting Eq. (11) into Eq. (4) yields the dynamical equations for the 
mode amplitudes 

fk,nn'{t) + J2(^\^kii)\^)fk,nn'{t) = 0, (14a) 
n 

where 



{n\ujl{t)\n) = J dh 4>*^ujl{h,t)4)n 



(14b) 



For our geometrical conditions and with the assumption that the initial 
charge configuration in Eq. (2) is the adiabatic vacuum, annihilated by ak,n{0) 
and by &fc,n(0)) the expectation value of the current becomes 

-{k-eA{b,t))J2Mb)rn'ib) 

(15) 

X Yl^fln'n{t)fkMt)+fk,n'flit)fk,nfl{t)]- 
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As is well known (see, e.g., the broad discussion in Ref. [17]), the current 
in Eq. (15) suffers from divergences. For the case of infinite volume, this has been 
dealt with systematically through the procedure of adiabatic regularization [17]. 
Unfortunately, that is inapplicable here because of the mode mixing inherent in 
our problem and explicit in Eqs. (14). We therefore follow an alternative strategy 
here. First we consider the high-momentum limit {k — > oo, n — > oo) for fk,nn'{t) 
satisfying the dynamical equations (14) and the wronskian condition (13). A direct 
generalization of the usual WKB approach [17] suggests 

fk,n'n{t) > {n I ^2iOk{t) 

as k, n, n' — > oo, where we restrict ourselves to the leading term of the WKB since 
the higher-order refinements are of no use in the presence of mode mixing. It is 
easily seen that this form satisfies Eq. (13) exactly and Eqs. (14) through order 
k^. To study the worst divergence in Eq. (15), we substitute this into the factor 
involving a summation over n, n', and n there to obtain, using closure, 

/•oo 1 

J2Mb)K'{b) ■ 27r / h'dh'<j>l{h')—-—ci>n'{h') = 

•^0 ^^^^b^t) ^^^^ 

where we again used closure to obtain the last line. It is then easy to see that this 
divergent form in fact vanishes upon symmetric integration [10,17] in Eq. (15). 
This identifies the worst, cubic divergence there, which is dealt with by subtracting 
the form of (17) to arrive at 



-{k-eA{b,t))Y,<Pnm*n'{b) 

^ 1^ '^^[fk,n'nfk,nn + fk,n'nfk,nn\ ~ cOk{b t) 



Ort.T}.' 



This expression still contains, of course, the usual logarithmic divergence 
of charge renormalization, dealt with in adiabatic regularization [10,17] by con- 
sidering higher-order WKB than shown in Eq. (16). This is what does not go 
through in the mode-mixed case. Instead, we base ourselves on the expectation 
that charge renormalization is what must emerge in any event. This is shown 
explicitly for adiabatic regularization in the infinite volume in Ref. [17], and we 
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do not expect a change in high-momentum, short-interval behavior here merely 
because of the presence of the cylinder walls. Thus we calculate Eq. (18) with a 
A;-momentum cutoff A and a radial mode cutoff Umaxj the same time replacing 
the charge e by the renormalized charge 



-1/2 

, (19) 



and testing that results are insensitive to A and Umax- 
Equation (18) is now used in Eq. (2) as 

^-viy{b,t)={jib,t)) (20) 

for the dynamics of the electromagnetic field. The coupled equations that must 
be solved for the field-theory study of back-reaction in an infinitely-long cylinder 
are then Eqs. (14), (18), and (20). 



III. TRANSPORT FORMALISM 

The transport equation we must deal with in the presence of the cylindrical 
symmetry of Figure 1 is, for a positive charge e, 

|/(b;p;<) + ^.^/(b;p;<) ^^^^ 

+ e (e + I X b) . (^ + S A)/(b; p; = SIE. B; p}, 

where / is the transport function and 5" is a source function based on the Schwinger 
mechanism for pair production in a fixed field. The spatial variable b is, as hith- 
erto, the transverse position measured from the cylinder axis; again there is no 
2-dependence here. The momentum p = {p_L, py} is taken in terms of a transverse 
component pj_ and a longitudinal one while £ = (p_L^ + + m^)^/^ is the 
particle energy. Note that with our geometry there can be dependence only on 
an azimuthal angle ip, the angle between b and p^. In parallel to the field case of 
the previous section, we impose the boundary condition f{h = R; p_L,p\\; t) = 0. 
Given the field configuration discussed in the previous section and shown in Figure 
1, the transport equation (21) takes on the form 



dp\\ \ £ J P\\ \^ J dp 



f = S. (22) 
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Negative charges move in accordance with Eq. (21) or (22), but with e replaced 
by — e. Since the source term depends only on |e|, it follows from Eq. (22) that the 
transport function for negative charges is given by that for positive charges with 
P\\ ^ -P\\- 

The source term in Eq. (21) or (22) is taken as [6-10] 

S{E, B- p} = 5{p\\)\eWE^-BHog 

where we insist that the invariant E"^ — B'^ be timelike, so that energy consid- 
erations allow pairs to be produced, taking 5' = for S > In Eq. (23), the 
usual choice [6-10,12] has been made to the effect that the particles are produced 
at zero longitudinal momentum. A study based on the projection method [13] 
has found that the source term is well localized around ~ for rather large 
fields, eVE'^ — B'^/^/p±'^ + m? > 1 in our notation; this situation is, in fact, the 
only one considered here. Since our results depend ultimately on the pair cur- 
rent, which involves an integration over py, they prove to be very insensitive to 
the appearance of 5{p\\) in S, being essentially unchanged for distributions in p\\ 
centered around p\\ = 0, normalized with respect to integration over as is S(p\\), 
and with widths ranging up to several times m. The nonmarkovian features found 
in Ref. [13] also disappear, of course, with the source term of Eq. (23), in con- 
sonance with the findings of Ref. [13] for large fields. Last, we shall not consider 
here the consequences of Bose enhancement in the transport equation treatment, 
the general eff'ects being well known from previous work [6,10]; to deal with such 
enhancement for the finite-volume case would require some selection of cell size in 
phase space within which statistical effects would take place. 

The nontrivial Maxwell equations in our geometry are 

dB dE BE .Id,, , , 

where j is the current — again purely longitudinal — given by a conduction and a 
polarization part, 

j {b, t) = jcond {b, t) + jpoi {b,t), (25a) 

where 

icond(6, t)^2ej ^^/(b, p; t) (25b) 
9 



1 + exp 



\eWE^-B^ 



(23) 



and 



(25c) 



The coupled equations that must be solved for the transport-theory study of back- 
reaction are then Eqs. (22), (24), and (25). 



We discuss first the procedures for the field-theory formulation. Since the 
three-dimensional studies of back-reaction in an infinite volume proved taxing from 
the computational point of view [10], it was, of course, to be expected that the 
finite- volume case would present difficulties in terms of the length of processor time 
required, and this was indeed the case. The infinite-volume calculation required 
very refined grids for the time variable (typically the time step was on the order of 
10~* m~^) and longitudinal momentum (typically requiring some 10^ points). The 
renormalization scheme used here [see Eq. (19) and the discussion surrounding it] 
is, however, simpler than the iterative one used for infinite volume [10]. 

Niimerical procedures were straightforward: the current j(6, t) and poten- 
tial A{h, t) and field E{h, t) were all expanded in terms of the basis set of Eq. (6). 
The initial field was usually taken in the form 



in the notation of Section II. The mode amplitudes fk,nn'{t) were then calculated 
from the coupled ordinary differential equations (14) using Runge-Kutta and the 
Fourier-Bessel components of the current j{b,t) were evaluated by Simpson inte- 
gration of Eq. (18). The subtraction of the counter term in (18) required a return 
to configuration space for the evaluation of a;j^^(6, t). The Fourier-Bessel compo- 
nents of the potential and field were then advanced using Eq. (20) — again an 
ordinary differential equation in the transform space — with Runge-Kutta. Thus 
worries about possible instabilities in partial differential equations were avoided. 

The rapid fluctuations encountered in the mode amplitudes in the inflnite- 
volume case seem here to translate into a great deal of mode coupling in the flnite 
volume, with a consequent need for a large number of Fourier-Bessel components 
in order to represent the current adequately. We here worked typically with about 
500 longitudinal momentum points and summed n out to about fimax ~ 

20. Time 



IV. NUMERICAL PROCEDURES 



E{b,0)^EQjQ{ZQQb/R), 



(26) 
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steps again had to be kept small {dt ~ 10~^). Since computation times rise 
quartically, or faster, with rimaxi this led to runs of several weeks with a 100-MHz 
or 150-MHz processor. Moreover, as for infinite volume, quantum fluctuations 
will eventually be more rapid than can be accommodated by any grid choice, 
so that our computation scheme must eventually break down. Confidence in the 
validity of the field calculation is therefore partly based on the usual tests checking 
that results do not change appreciably when grid sizes are refined. It also partly 
derives from the qualitative and even semiquantitative agreement between field 
and transport methods over the time range we treat, a comparison which is, of 
course, the primary goal of this study. 

Somewhat surprisingly, the transport formalism is also somewhat tedious 
computationally. The infinite-volume case allowed for solution of the transport 
equations using essentially an analytic scheme based on the method of character- 
istics [10], but for a finite volume we had to solve Eq. (22) as a partial difi^erential 
equation, using a Lax method [18]. Other methods, such as Lax-Wendroff or stag- 
gered leapfrog seemed to offer little improvement over this. The time evolution 
of the electromagnetic fields as partial differential equations (24) — even in the ab- 
sence of a current j — proved unstable and this was not easily cured through the use 
of Lax, Lax-Wendroff, or leapfrog methods, so we again resorted to Fourier-Bessel 
decomposition and advanced the coefficients in time using a predictor-corrector 
method. The time increment was again dt ~ 10~^ and roughly 120 longitudinal 
momentum points were required. In the transverse direction about 120 momentum 
points and 80 spatial points were needed. 

The outcome of all this was again very long computation time. The situ- 
ation was eased somewhat by great insensitivity to the treatment of the angular 
variable Eventually we chose to average over this, fixing (/? = 7r/2. This changed 
the results typically by only a few percent. We also eventually averaged the trans- 
verse momentum |pj_| in a manner suggested by the gaussian dependence of the 
source term of Eq. (23), taking (p^) = {^n\e\/2){E'^^B'^Y^'^. This approximation 
reproduced very well the amplitudes and frequencies of the oscillatory behavior 
of j{b,t) and E{b,t), but was clearly deficient near the zeros of these quantities, 
where it led to a flattening or step in j near its zeros. 
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V. RESULTS AND DISCUSSION 



In presenting results in Figures 2 through 13, all quantities are given in 
natural units, so that spatial coordinates are measured in 1/m and momentum 
coordinates in m. The electromagnetic four-potential is presented as a{b, t) = 
eA{h, t)/m^, the electric field as e(6, t) = eE{h, t)/rn^, and the current as j{b, t) — > 
ej{b,t)/m^ . The cases chosen for presentation here parallel parameter sets picked 
in the study of the infinite- volume situation [10]. We shall generally show the 
current and the electric field averaged over the cylinder cross section according to 



where q{h, t) = j{h, t) or e(b, t). This, as we shall see, converges much faster as the 
number of Fourier-Bessel components rimax is increased than do, say, the on-axis 
values. 

Figures 2 and 3 present the transversely-integrated current and electric 
field for unrenormalized charge = 10, on-axis initial electric field e{b = 0, t = 
0) = 6, and cylinder radii R = 2 and 5. The large value of the charge chosen 
here has the advantage, given our need for very lengthy computer time, that back- 
reaction effects have a chance to make themselves felt within the relatively short 
time interval t < 4 (sec Ref. [10]). As will be generally true below, the electric 
field integrated over the cylinder cross section is very similar for the field-theory 
approach and the transport calculation. Of course, this is partly the case because, 
for the small radii considered here, the electric field oscillations are largely driven 
by cylinder radius, not by back-reaction. Thus, by i? = 10 this agreement begins 
to break down. On the other hand, R = 10 stretches the range of reliability of 
our calculation, since we must require that the "transverse momentum" rimaxT^ /R 
be on the order of 10 to 20 in order to cover the relevant momentum distributions 
and to yield a valid regularization scheme using Eq. (18). Since in practice we 
are limited to rimax < 20, this implies a breakdown in our procedures for R much 
beyond 5 or so. The transversely-integrated currents as calculated from field theory 
are substantially more oscillatory than their transport-theory counterparts, which 
is not surprising in view of the infinite- volume results [10]. The initial sharp rise 
in the current occurs in both cases, though with rather different values; thereafter 
there is reasonable similarity between the two methods if one takes into account 
an effective time-averaging over the field-theory fluctuations. It is worth noting 




(27) 
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[19] that the initial rapid variations in the electric field in the quantal case may 
"shake off" charges by means of the fast fluctuations, an effect not included in the 
transport approach here, where only pair tunneling is included. 

In an attempt to make some connection with the infinite- volume results, 
we show j and e on the cylinder axis for = 10, initial field e{b = 0, t = 0) = 6, 
and R = 2 and 5 in Figures 4 and 5. Agreement here between field theory and 
transport formalism is much less satisfactory than the situation for the integrated 
values: The electric field on-axis remains similar for the two cases, again because its 
values are essentially determined by the initial condition and the cylinder radius, 
but the current differs by nearly two orders of magnitude. This is almost certainly 
in large part a numerical breakdown brought about by the difficulty in using a 
large enough number of Fourier-Bessel components to achieve reliable results for 
the on-axis case. (For Figures 2 through 5 we used Umax = 20.) Since Jo(0) = 1, 
all the components add coherently to produce the 6 = values, and thus these are 
very sensitive to limitations in Umax- On the other hand, the on-axis value of the 
current is weighted with a factor of b and hence does not really contribute when 
global features are considered. 

The corresponding picture for a smaller charge, namely, = 4, with ini- 
tial field e{b = 0, t = 0) = 7, and again for R = 2 and 5 is shown in Figures 6 
through 9 (where we have worked with Umax = 10 in order to allow for the calcu- 
lation to extend out to t = 30, and encompass several cycle times). As is known 
from the infinite- volume case [10], the field-theory and transport-formalism results 
eventually drift out of phase, a condition that would be improved by incorporating 
Bose enhancement into the transport equation [10]. Again, if one were to average 
over the fiuctuations present in each cycle of the current, agreement between the 
two methods is quite reasonable for the integrated values, and not very good for 
on-axis results. 

In Figures 10 and 11 we compare, for = 10, e{b = 0, t = 0) = 6, 
and R = 2 and 5, results for two different numbers of Fourier-Bessel compo- 
nents, Umax — 10, dt = 2 X 10~^, 401 fc-grid points, and A = 20 (dashed curve) 
vs. Umax — 20, dt = lO""^, 1,001 /c-grid points, and A = 40 (solid curve) for 
the transversely-integrated current and electric field. The same comparisons for 
the on-axis quantities are shown in Figures 12 and 13. Since we have here a 
doubling in this upper limit for the number of transverse modes rimax, the rea- 
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sonable agreement for the integrated values suggests that the results there for 
T^max = 20 are quantitatively reliable. Obviously this is not yet the case for the 
on-axis quantities — notably the current, which differs by two orders of magnitude 
for t > 2 — which at best have qualitative validity. It is perhaps not surprising that 
the less numerically reliable field-theory result (dashed curve) more closely resem- 
bles the transport calculation, since its inferior numerics automatically perform 
some sort of average over the quantal fluctuations of the field-theory result. Such 
an averaging would otherwise only properly be performed by doing the field-theory 
computation with much more refined parameters than are presently accessible and 
then averaging these results over time systematically, or — better still — averaging 
over momentum distributions fk,nn'{t) before going on to compute final physical 
values [10]. 

In sum, the present results indicate that back-reaction in a finite volume 
allows for the replacement of the full field-theory calculation by the much simpler 
transport formalism at least at a loose, qualitative level for quantities averaged over 
the transverse direction. It is perhaps not surprising that the very close linkage 
between field theory and transport formalism found in the case of infinite volume 
[6-10] is lost here: The values of the radius R that we deal with are comparable 
with the particle mass, so that quantal effects should be expected to enter, and the 
spatial variability of the field is such that the "local" Schwinger form of Eq. (23) 
is not likely to prove adequate. This is in contrast to the infinite-volume situation 
[6-10], where the temporal variability of the physical quantities is measured by 
cycle times roughly an order of magnitude larger than R here. 

This rough, qualitative correspondence between field theory and transport 
formalism may also apply for the physical quantities in greater detail as functions 
of the transverse parameter 6, but present computer capabilities do not yet allow 
one to establish this. In all likelihood the correspondence will be strengthened 
somewhat by the consideration of the entire physical context, including particle 
interaction, irregular geometry, and the like. Thus for many practical purposes one 
may be able to replace a field-theory approach by a transport one at the level of 
qualitative behavior, and use the latter for further study of pre-equilibrium parton 
production. 
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Figure captions:- 



FIG 1. Cylindrical geometry considered here. On the left side of the cylinder are 
shown typical actions of forces on negative and positive charges. Since charges 
of opposite sign are evenly distributed and their transverse motion is the same, 
they neutralize each other insofar as transverse contributions to the current are 
concerned. 

FIG 2. Comparison between field-theory results (left-hand side) and transport 
theory calculation (right-hand side) for quantities integrated over the cylindrical 
cross section [i.e., tranversely — see Eq. (27)]. Here the charge is = 10, the 
initial electric field on the cylinder axis is E{b = 0, t — 0) — 6, and the cylinder 
radius is i? = 2, all in natural units (see Section IV). 

FIG 3. Same as Figure 2, but for a cylinder radius R = 5. 

FIG 4. Comparison between field-theory results (left-hand side) and transport 
theory calculation (right-hand side) for quantities on the cylinder axis. The charge 
is = 10, the initial electric field on axis is E{b = 0, t = 0) = 6, and the cylinder 
radius is i? = 2 in natural units (see Section IV). 

FIG 5. Same as Figure 4, but for a cylinder radius R = 5. 

FIG 6. Comparison between field-theory results (left-hand side) and transport 
theory calculation (right-hand side) for quantities integrated over the cylindrical 
cross section. Here the charge is = 4, the initial electric field on the cylinder 
axis is E{b = 0, t = 0) = 7, and the cylinder radius is i? = 2 in natural units. 

FIG 7. Same as Figure 6, but for a cylinder radius R — 5. 

FIG 8. Comparison between field-theory results (left-hand side) and transport 
theory calculation (right-hand side) for quantities on the cylinder axis. The charge 
is = 4, the initial electric field on axis is E{b = 0, t = 0) = 7, and the cylinder 
radius is = 2 in natural units (see Section IV). 

FIG 9. Same as Figure 8, but for a cylinder radius R = 5. 

FIG 10. Comparison between computation for integrated J{t), A{t), E{t) [see 
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Eq. (27)] using time interval dt = 10~^, 1,001 longitudinal grid points, 20 Fourier- 
Bessel components {rimax = 20), and cutoff A = 40 taking several weeks (solid line) 
with a less accurate version having dt = 2x 10~^, 401 grid points, rimax = 10, and 
A = 20, taking only a few days. The physical parameters here are the same as for 
Figure 2. 

FIG 11. Same comparison as for Figure 10, but with cylinder radius R = 5. 

FIG 12. Same comparison as in Figure 10, but for on-axis quantities j{b = 
0, t), a{b = 0, t), e{b = 0, t). 

FIG 13. Same comparison as for Figure 12, but with cylinder radius R — 5. 
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